Optimization of alloy-analogy-based approaches 
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An analytical expression for the self-energy of the infinite-dimensional Hubbard model is 
proposed that interpolates between different exactly solvable limits. We profit by the combi- 
nation of two recent approaches that are based on the alloy-analogy (Hubbard-III) solution: 
The modified alloy-analogy (MAA) which focuses on the strong-coupling regime, and the 
Edwards-Hertz approach (EHA) which correctly recovers the weak-coupling regime. In- 
vestigating the high-energy expansion of the EHA self-energy, it turns out that the EHA 
reproduces the first three exactly known moments of the spectral density only. This may be 
insufficient for the investigation of spontaneous magnetism. The analysis of the high-energy 
behavior of the CPA self-consistency equation allows for a new interpretation of the MAA: 
The MAA is the only (two-component) alloy-analogy that correctly takes into account the 
first four moments of the spectral density. For small U, however, the MAA does not repro- 
duce Fermi-liquid properties. The defects of the MAA as well as of the EHA are avoided 
in the new approach. We discuss the prospects of the theory and present numerical results 
in comparison with essentially exact quantum Monte Carlo data. The correct high-energy 
behavior of the self-energy is proved to be a decisive ingredient for a reliable description of 
spontaneous magnetism. 
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I. INTRODUCTION 

A central problem in solid-state physics concerns in- 
teracting electrons on a lattice. In dealing with itin- 
erant magnetism, heavy-fermion compounds or high- 
temperature superconductivity, for example, electron- 
correlation effects are of exceptional significance. Much 
insight into the fundamental role of electron correlations 
can be gained by studying the Hubbard model |l| |[ || . 
In spite of its apparent simplicity, an exact solution for 
the whole parameter range is not available up to now, 
and a completely satisfactory understanding of its prop- 
erties has not yet been achieved. 

A notable exception is the extreme case of the one- 
dimensional model. In particular, the exact solution for 
the ground state is known in d = 1 B. The opposite 
limit of high spatial dimensions d, which has been intro- 
duced by Metzner and Vollhardt B, is likewise impor- 



tant. For d = oo there are considerable simplifications 
that are due to the momentum independence of the elec- 
tronic self-energy 0. However, the infinite-dimensional 
model still remains non-trivial. It is of special interest 
since its essential properties are expected to be compa- 
rable to those at low dimensions d — 2, 3. If an ap- 
proximation scheme was available that is reliable for the 
entire range of the model parameters in d — oo, this 
would provide a proper dynamical mean-field theory in 
any dimension d > 1 fjl @|. 

Quantum Monte Carlo (QMC) §, [To[ [n| and exact 
diagonalization methods (ED) |l2|, 13| can yield essen- 
tially exact results but also suffer from limitations: ED 
is restricted to a rather small number of orbitals, and 
thus a smooth density of states cannot be obtained. On 
the other hand, QMC yields its results for the discrete 
Matsubara energies or along the imaginary time axis. 
Therefore, it is difficult to access the low-temperature 
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regime where statistical errors become important within 
the QMC method. Furthermore, to obtain dynamical 
quantities it becomes necessary to perform an analyti- 
cal continuation to the real axis which constitutes an ill- 
conditioned numerical problem. For these reasons the 
development of analytical (but approximate) methods 
for the infinite-dimensional Hubbard model still remains 
necessary. 

An analytical approach to the d — oo Hubbard model 
should be guided by rigorously solvable limiting cases. 
There are a number of exact results which impose strong 
necessary conditions on any approximation: (i) Second- 
order perturbation theory around the Hartree-Fock so- 
lution @ (SOPT-HF) or self-consistent SOPT [§, |l| 
yield the correct self-energy in the weak-coupling regime, 
(ii) The exact low-energy (Fermi-liquid) properties of the 
self-energy and the Luttinger sum rule |[6[ in particular 
have been analyzed in Ref. || for the d — oo case, (iii) 
The atomic limit of vanishing hopping t = represents 
an extreme but important limit. The first non-trivial 
correction to the atomic limit can be obtained from a 
canonical transformation of the Hubbard model that has 
first been considered by Harris and Lange ji"?], For 
U i — > oo the approach provides rigorous information on 
the centers of gravity as well as on the spectral weights 
of the lower and upper Hubbard bands, (iv) The high- 
energy behavior of the self-energy is determined by the 
moments of the spectral density which can be calculated 
independently. For any U the high-energy behavior is 
important for a qualitatively correct description of the 
charge-excitation peaks. For U i— > oo the first four mo- 
ments are necessary to be consistent with the results of 
Harris and Lange (v) Finally, the Falicov-Kimball 
model (FKM) |2^] can be considered as a special limit 
of the Hubbard model. It is obtained by formally intro- 
ducing a spin-dependent hopping and allowing only one 
of the two spin species to be mobile. The d — oo FKM is 
exactly solvable as has first been shown by Brandt and 
Mielsch |2l[] . The self-energy is given by the Hubbard-III 
alloy-analogy (AA) solution p2| . 

It may be interesting to have a theory at one's 
disposal that for arbitrary fillings recovers the weak- 
coupling, the atomic and the FKM limit, that is able to 
reproduce the first four moments of the spectral density 
and thereby the exact strong-coupling results of Harris 
and Lange, and that allows to directly treat zero tem- 
perature and real energies. The construction of such an 



approach is the main purpose of the present paper. To 
make contact with the FKM, the approach will be based 
on the alloy-analogy solution. The AA represents an at- 
tractive approximation which is simple, self-consistent 
and non-perturbative. On the other hand, there are de- 
fects that show up in other limits and that have given 
rise to modifications of the conventional AA: 

A mean-field theory has been constructed in Refs. 
p^ , pi[ which is conceptually similar to the AA but in 
contrast to the latter derivable from an explicit free- 
energy functional and thus thermodynamically consis- 
tent and of variational character. This approach as 
well as the A A itself, however, cannot reproduce the 
aforementioned weak- and strong-coupling limits. The 
weak-coupling behavior of the AA is corrected within 
the Edwards-Hertz approach (EHA) j2|, f|| §7], |§]. A 
drawback of the EHA and also of the AA, however, con- 
sists in the fact that the strong-coupling results of Har- 
ris and Lange cannot be recovered. A recently proposed 
modified alloy analogy (MAA) improves upon the 
AA and ensures consistency with the results of Harris 
and Lange for U i— ► oo. On the other hand, the MAA 
is not able to reproduce the Fermi-liquid properties for 
small U . 

Our study aims at a proper combination of the 
MAA and the EHA which keeps their advantages and 
avoids their defects. After the necessary theoretical 
preparations (Sec. II), we derive a helpful alternative 
interpretation of the MAA in Sec. III. Sec. IV concerns 
the high-energy expansion of the EHA self-energy. 
The analysis of the MAA and the EHA is a necessary 
condition for the concise presentation of the combined 
theory in Sec. V. Numerical results obtained within 
the new method will be discussed for the paramagetic 
phase and in particular with respect to ferromagnetism 
(Sec. VI) before we come to the conclusions in Sec. VII. 

II. PREPARATIONS 

Using standard notations, the Hubbard Hamiltonian 
reads: 

H = {Tij - uSij) c\ a c j(T + ~U^2 n irT m- a . (1) 

ija ia 

The limit d i— > oo has to be taken with the proper scal- 
ing of the hopping integral between nearest neighbors to 
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ensure the model to remain non-trivial || . For a hyper- 
cubic lattice with coordination number Z — 2d and hop- 
ping T(jj) =t = t* I \/2Z (t* — const.) we get a Gaussian 
for the (non-interacting) Bloch-density of states (BDOS) 



1 



-(E/t*) 2 



t* 



(2) 



A generalization of the d = 3 fee lattice to infinite dimen- 
sions with Z = 2d(d— 1) results in a strongly asymmetric 



Bloch-density of states 1 30 



e -(l + v / 2S/t*)/2 

t*Jn{l + V2E/t* 



(3) 



where the hopping integral between nearest neighbors 
T tj is positive and t = t*/VZ 

In infinite dimensions the self-energy Ti a (E) is k 
independent or site-diagonal ||. For a homogeneous 
phase, the on-site one-electron Green function thus de- 
pends on the lattice geometry via the BDOS only: 



J —< 



E - (z - fi) - 



dz . 



(4) 



Let us consider its high-energy expansion: 

^(*0 = »£^T- (5) 



m— 



The expansion coefficient are the moments 
1 f°° 

M^ = - E m A a (E)dE (6) 

J — oo 

of the spectral density A a (E) = -~lvaG a (E + i0+). 
They can be calculated exactly by making use of an 
alternative but equivalent representation: 



M M H[£ „ w t J+) 



(7) 



Here CO — [O, H]- denotes the commutator of an op- 
erator O with the Hamiltonian, and [•••,•■■]+ is the an- 
ticommutator. The calculation along Eq. (]?]) is straight- 
forward. However, with increasing order m the moments 
include equal-time correlation functions of higher and 
higher order which have to be expressed by some means 
in terms of known quantities again. This fact limits the 
number of moments that can be used in practice. For 



the Hubbard model the moments are useful up to m = 3 
(see Refs. [|32[ p3|, for example). In particular, we obtain 
for d = oo: 



M< 3 > 



= 1 



T + U ti — (j , 

^TijTji + 2T Un_ a + XJ 2 n_ a , 



jk 



TijT jk T ki + 3U Th — ^ ^ TijTj. 



+ T U 2 n- a (2 + n-a) 
+ U 2 n- a {l-n- a )B- 



U 3 n- a 



For abbreviation we have defined: = 



(8) 

-/j,6ij, To — 

Tu and B a = B a — //. This so-called 'band-shift' consists 
of higher-order correlation functions, 

1 



B„ 



T 



£^■(4^(2^-1)), ( 9 ) 



but can be expressed in terms of the one-particle Green 
function |33|]. Note that while the hopping Ty as well 
as the correlation functions (c\ cr Cj (7 (2n i _ (7 — 1)) scale as 
1/y/d for d t-y oo, the lattice sum in (§) remains finite. 
Contrary to the finite-dimensional case, however, the 
band-shift is k independent for d = oo. 

The coefficients in the 1/E expansion of the self- 
energy, 



E m 



(10) 



are obtained by inserting Eqs. (|^) and (]Tc ) into (Q) and 
taking the moments of the BDOS from (|) for U = 0: 

C(°) = Un-„ , 

= U 2 n_ a (1 - 7i — fj ) , 

= U 2 n^ (1 - n_ ff ) (B_ a + U(l - n_ ff ) 



(11) 

Let us mention some available rigorous results for 
the weak- and the strong-coupling regime. For U i— > oo 
the spectrum is dominated by the two charge-excitation 
peaks (Hubbard bands). At each k point the centers 
of gravity as well as the spectral weights of the lower 
and the upper Hubbard band can be calculated ex- 
actly within a perturbational approach due to Harris 
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and Lange |17], Q that uses t/U as an expansion pa- 
rameter. For d = oo we have: 

T^ L) (k) = (1 - n_ CT )e(k) + Ti—tjB—a + 0{t/U) , 
T 2 ( ^ L) (k) = U + n^e(k) + (1 - n^)B. 
a£ L) (k) = l-«r'W = l-«-« 



■WO, 

0(t/CT) . (12) 



The only difference to the finite-dimensional case Il7j 
consists in the fact that the k dependence of the centers 
of gravity Tp^ L \h) and the weights a-p^ (k) (p = 1, 2) 
is exclusively due to the Bloch dispersion e(k). For finite 
d an additional k dependence is introduced via the band 
shift. 

In the weak-coupling regime the usual perturbative 
approach applies. Up to order U 2 the self-energy is given 
by: 

4 SOPT > (E) = Un_ a + 4 soc > (E) . (13) 

The first-order term is the Hartree-Fock self-energy; the 
second-order contribution reads: 

(f(x)f(~y)f(z) + f(-x)f(y)f(-z)) dxdydz(U) 

Here f(x) is the Fermi function, and A^(E) is defined 
as the free (U = 0) spectral density shifted in energy by 
a (possibly spin-dependent) constant E a : 



A^{E)=A^\E-E a ) 



(15) 



The shifts E a are introduced for later use. With E a = 
the plain or conventional second-order perturbation 
theory (SOPT) and with E a = Un^ a the SOPT around 
the Hartree-Fock solution E3 is recovered. Both 
versions of SOPT as well as the self-consistent SOPT 
M are identical up to order U 2 . 



III. MODIFIED ALLOY- ANALOGY 

The original alloy-analogy solution (AA) by Hubbard 
p2| is at variance with both, the exact weak- and strong- 
coupling results. Recently, a modified alloy- analogy so- 
lution (MAA) has been proposed (2^] to overcome the 
drawbacks of the conventional AA with respect to the 
limit XI i — ► oo. Any two-component alloy analogy re- 
quires the specification of the two atomic levels and the 



corresponding concentrations. The self-energy is then 
obtained from the CPA equation [p4[: 



E pa — S CT (E) — T 



-G a {E) [E pa -V a (E)-T ] 



(16) 



Within the conventional AA the levels and weights are 



E 



(AA) 



Ict 



taken from the atomic-limit solution J 
T , E^ AA) =T q + U and x[ AA) = 1 - , x£ A) = ri- _ 
This choice, however, is by no means predetermined. 
Within the MAA [§§ the self-energy E^ MAA) (E) is still 
obtained from the CPA equation (|l6|), but the atomic 
levels and the concentrations are replaced by: 

4^ AA) = \ (To + U + 5_ ff ) + (-If x 



i (U + B_ a - T ) 2 + Un_ a (T - B_ a ) , 



and 



(MAA) 



U(l 



E 



(MAA) 



„(MAA) pfMAA) 



1 



(17) 



(MAA) 
l 2er 



(18) 

The expressions include the band shift B- a . If B- a 
is replaced by T , the MAA reduces to the conventional 
AA. As can be seen from Eq. (^|) this replacement is cor- 
rect for the atomic and for the "FKM limit". Further- 
more, -B-jj = To is required by particle-hole symmetry 
in the case of a paramagnet at half-filling. In all other 
cases the MAA is different from the AA. In particular, 
via the band shift the atomic levels might now become 
spin-dependent. 



In the original work 1 29 the MAA has been derived 
by referring to the spectral-density approach (SDA) 
p^ , |33| . p5{ . Considering the so-called split-band regime 
of the CPA Q, the consistency of the MAA with the 
results of Harris and Lange (|l^) could be proven for 
U h-> oo. 

For later purposes it is important to give another 
interpretation of the MAA which is different from the 
original one in Ref. |2{J: Since the atomic values and 
the concentrations that are needed to construct an alloy 
analogy for the Hubbard model are not predetermined, 
one can ask for the optimal choice (in some sense). In 
the following we argue that the MAA is obtained if one 
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demands that the resulting theory should correctly take 
into account the first four moments of the spectral den- 
sity given in Eq. (|J): 

We consider the CPA equation ( |l6| ) as an equation 
for the self-energy that contains four unknown param- 
eters, E p(J and x pa (p = 1,2), which have to be fixed 
by imposing the correct results for the moments of the 
spectral density. Since the moments are closely related 
to the high-energy behavior of the Green function, wc 
have to expand the CPA equation in powers of 1/E us- 
ing Eqs. Q and (10). Considering terms up to l/E 2 , 
this yields the following equations: 



J2- 



'pa , 



p 

[(^-T -C(°)) 2 Af(°)-CW 

p 

$> pff [{E pa -T,-C^f{M^f 
p 

+ (E pa -T -CWfMW 

- 2(E pa -T - Cj?>)CWMW - C«l . (19) 



Inserting the exact expansion coefficients (|8|) and (|ll|) for 
the Green function and for the self-energy and solving 
for Ep a and x pcr , indeed yields the MAA results (0) and 
©■ 

This proves that the MAA not only recovers the ex- 
act strong-coupling results of Harris and Lange, but also 
(for arbitrary U) respects the conditions that are im- 
posed on the overall shape of the spectral density by 
their exactly known first four moments. In this sense 
the MAA can be termed to be an optimal alloy analogy. 

It is not very difficult to analyze the high-energy be- 
havior of the conventional AA. Starting from Eq. ( |lS| ) 
and inserting the AA levels and concentrations, it shows 
up that the last equation cannot be satisfied. While 
the AA yields the correct coefficients up to order 1/E 
in the expansion of the self-energy, it turns out that 

(2) 

the coefficient is not recovered. Instead, one gets: 
Ci 2 ,AA ' = |b_ <t h^t - This is only exact in the atomic 
limit. Equivalently, it can be stated that the AA re- 
covers three (instead of four) moments of the spectral 
density only. 



While the MAA (in the sense of Harris and Lange) 
is correct in the strong-correlation regime, there is a 
severe drawback of the method: It fails to reproduce 
the Fermi-liquid properties for small interactions U. 
The same defect is inherent in the conventional AA. 



IV. EDWARDS-HERTZ APPROXIMATION 

A correction of the conventional AA with respect 
to Luttinger's requirements and Fermi-liquid properties 
is due to Edwards and Hertz p5| . Just as the MAA 
improves upon the strong-coupling regime of the AA, 
the Edwards-Hertz approximation (EHA) corrects the 
A A in the weak-coupling regime. 

The self-energy within the Edwards-Hertz approxi- 
mation (EHA) is implicitly given by p5l (see also Refs. 

MM)-- 



i(EHA) 



(E) = 



Un- a 



l-\G a {E) (17 -E 



(EHA) 



(20) 



This equation is formally equivalent to the CPA equa- 
tion ( |ilf ) provided that the atomic levels and the concen- 
trations are taken to be those of the conventional AA. 
Contrary to the AA, however, the one-particle Green 
function that occurs in the CPA equation, G a (E), has 
been replaced by: 



G a (E) 



t/ 2 n_ CT (l 



Wsoc) (£; _ e ^eha) ( £ )+ £ ct) 

-a) 

(21) 

£,(SOC) j g second-order contribution to the SOPT 
self-energy introduced in Eq. (|l4|). The replacement is 
more or less ad hoc to enforce the correct weak-coupling 
behavior of the theory. 

The shifts E a that appear in the definition ( |l5| ) of 
the spectral density (E) and once more in Eq. (|2l| ) 
have been introduced by Wermbter and Czycholl [ p7j . 
They are necessary to ensure that the EHA is exact in 
the atomic limit for all band-fillings. The shifts have to 
be determined self-consistently by the following require- 
ment: 

/oo />oo 
f(E)A a (E)dE= / f{E)A^{E)dE. 
-oo J — OO 

(22) 

The EHA should be regarded as an interpolation 
scheme. It yields the correct results in the atomic as 
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well as in the FKM limit and is exact for small U up 
to order U 2 |p7| . The latter implies that the expansion 
in U satisfies Luttinger's requirements. However, Fermi- 
liquid behavior breaks down in a parameter region where 
the series diverges p5[ . 

Let us now investigate the high-energy behavior of 
the EHA self-energy: 



4 EHA) (£) = E 



oo „(m,EHA) 



m— 



(23) 



First, we need the high-energy expansion of the modi- 
fied propagator G a {E). Applying the Kramers-Kronig 
relation, it can be written in the form: 



G a (E) = 



U 2 n- a (l - n_ CT ) 

-ilmSi SOC) (£;' + iO+) 
E-Y, i ® HA) (E) + E (T - E' 



dE' . (24) 



When expanding in powers of l/E, the following rela- 
tions are needed: 

J -hmJ:f oc) (E + iO + )dE = U 2 n- a (l - n-„) , 

J -hm SS S ° C) (E + *0 + ) E dE = 

U 2 n_ a (l-n_ a ) (B { V-n + E a 



(25) 



Using Eq. (|22|) these results can be derived from Eq. (jlj) 
immediately For abbreviation we have defined: 

BP = To + , 2 "7 g "\ E^<4^) (1) , (26) 



n a (l - n a ) 



where the one-particle correlation function results from 
the free (U = 0) but shifted off-diagonal Green function: 



{CLC 



h I f(E)lmG\jl(E + iO + ~ E a )dE. The 
expansion of the modified propagator then reads: 

\G a {E?) = 1 + (b^I - M + C/n_ CT ) + • • • . (27) 

The coefficients in this expansion are the moments of 
the (fictive) spectral density A a {E) = — ^Im G a (E + 
i0 + ) (This will become important in the next section). 



Inserting the result and the ansatz (|23|) in Eq. ( |20| ) and 
grouping the terms in powers of l/E, we finally obtain: 

C (0,EHA) = Un _ a 

(l-n_ a ), 



C (1,EHA) = t/ 2 n _ 



p(2.EHA) _ jj2 



B 



(i) 



r (1 - n_ CT ) x 
/i + U(l — n- a 



(28) 



Comparing with Eq. (^_1J), we notice that the 
first two coefficients ci°^ and Ca^ are predicted 
correctly. As in the conventional AA, however, the 
third expansion coefficient turns out to be wrong: 



(2, EHA) 



c, 



(2) 



Again, the replacement 



B^ a is correct in the atomic limit only. Fur- 



a 

B- 

thermore, via Eqs. (^), (||) and ( ^0[ ) we can conclude 
that the fourth moment of the spectral density, 
is not recovered within the EHA. Finally, this implies 
that for U i— > oo the EHA is at variance with the exact 
results of Harris and Lange [Eq. jl^)]. 

V. COMBINING THE APPROACHES 

We are now in a position to profit from the preced- 
ing analysis of the MAA and the EHA. We like to tackle 
the question whether a proper combination of both ap- 
proaches is able to avoid their respective defects and to 
keep their advantages. 

The following simple idea suggests itself: We start 
from the general CPA equation ([l6]), replace the one- 
particle Green function G a (E) by the modified propa- 
gator G a {E) as in the EHA, but take the atomic levels 
E pa and the weights x pa from the MAA according to 
Eqs. @ and @. It can be shown, however, that such 
an ansatz is neither correct up to order U 2 nor it recovers 
the high-energy coefficients up to order l/E 2 . 

We thus have to disregard the MAA levels and 
weights, but can try to retain the basic concept of the 
approach. In Sec. IV we have interpreted the MAA as 
being the only (two-component) alloy analogy that cor- 
rectly takes into account the first four moments of the 
spectral density: Insisting on the correctness of the mo- 
ments, unambiguously determines the levels and weights 
that appear in the CPA equation. Therefore, if we 
change the CPA equation by introducing the modified 
propagator G a (E) , we must allow the parameters to ad- 
just in order to keep the moments correct. This concept 
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turns out to be successful indeed. For clarity in the no- 
tations let us refer to this approach as an "interpolating 
alloy-analogy-based" approximation (IAA) since we can 
show that it recovers the weak-coupling limit and sat- 
isfies the requirements of Harris and Lange for strong 
U. 

To begin with, we once more consider Eq. ([l6]) with 
G a (E) replaced by G a {E) and with a priori unknown 
parameters E pa and x pa . As in the EHA the modi- 
fied propagator G a (E) is defined via Eq. ( pl| ) . To make 
contact with the moments of the spectral density, we 
have to expand the equation in powers of 1/E. This 
yields Eq. ( |l9| ) where the quantities now have to 

be interpreted as the moments of the modified propa- 
gator G a (E). Fortunately, only the first two moments 
are needed which can be taken from Eq. (|27]) . To ensure 
the correctness of the moments (|J) within the IAA, we 
then have to insert the exact expressions (^Tj) for the 
coefficients Ca m \ which yields the following system of 
equations: 



.(IAA) 



1 



$>£ AA) « AA) " To) = t/n_ CT , 
p 

£4 IAA > (E^ - T ) 2 = U 2 n^ 
v 

£4 IAA > « AA ) - To) 3 = C/ 3 n_ CT 
p 

+ C/ 2 n_ CT (l - - Si 1 ] ) . (29) 



Solving for E p 1 ^^ and Xp^^ , we obtain: 
EjW = T + \ (U + B_ a S«) + (-1) 



\ (U + B- a - B { H) 2 + X]n- a (si 1 ] - ) 



(30) 



and 



„(IAA) 



B- a - + To + U(l - 71 — fj ) - E\ 



(IAA) 



p(IAA) ^(IAA) 



= 1 -X 



(IAA) 



(31) 



The results turn out to represent a slight modification 
of the corresponding MAA results which are obtained if 
B^l is (ad hoc) replaced by T . In the last section we 
have recognized that the more simple one-particle corre- 
lation functions B^ a instead of the higher-order correla- 
tions -B_cr appear in the expansion coefficient c£ 2 ' EHA ^ 
of the EHA. Thus, it can be well understood that the 
difference B- a — B^l is contained in ( |30| ) and (|3|). If 
the difference is neglected, we will get the atomic levels 
and weights of the conventional AA again. 

This concludes the theory. However, we still have to 
check whether the IAA self-energy is correct up to order 
U 2 . This is by no means obvious since the U dependence 
of the atomic levels and the weights has become more 



complicated. Using the results ( |30| ) and (31) in the CPA 
equation (|l6|), we obtain for the IAA self-energy: 



s aAA) (jB) = 



j-G a (E) (U-Z£ AA) (E) 
1 - j-G a {E) (-B-cr - B 



(32) 



The expression clearly shows the differences with respect 
to the conventional AA as well as to the EHA. The equa- 
tion can be expanded in powers of U up to U 2 without 
knowing about the explicit U dependence of the modi- 
fied propagator. However, we need: 



B_ r 



B 



(i) 



const, x U + 0(U 2 ) , 



(33) 



which can be verified easily Using (p3[), it turns out 
that up to order U 2 we have: 

4 IAA )(£) = \J yl — fj + U 2 n_ a {l - n_ a )\G a {E) , (34) 

n 

which is the same result as found in the EHA and thus 
proves the IAA to be exact up to order U 2 indeed (cf. 
Ref. 0): 

S £AA) = Un _ a + s (soc) +0 {U 3 ). (35) 

Let us briefly check the other limiting cases men- 
tioned: It goes without saying that the trivial cases 
of U = 0,n = 0orn = 2as well as the requirements 
imposed by particle-hole symmetry and the Herglotz 
properties |3j| are fulfilled. Suppressing the hopping of 
the .[-electrons, we have B^ = B^ = T from Eqs. (||) 
and (Eq). Hence, the theory reduces to the EHA and 
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is thereby |27|] correct for the Falicov-Kimball model 
and also in the atomic limit. By construction, the 
IAA self-energy has the correct high-energy expansion 
coefficients up to order 1/E 2 , and thus the moments of 



the spectral density up to Al!f are reproduced 

exactly. It can be shown analytically (similar to the 
discussion in Ref. ]l9| ) or can simply be taken from 
the results presented in the next section that for strong 
U the k-resolved spectral density for each k-point 
consists of two dominant peaks (Hubbard bands). 
Their energetic distance is roughly given by U. Since 
the first four moments of the spectral density are also 
conserved for each k-point separately |3^] , the centers of 
gravity as well as the weights of both, the lower and the 
upper Hubbard band, are fixed and can be calculated 
analytically for U i— > oo along the lines of Ref. ]3^] , for 
example. The result turns out to be identical to what 
has been derived within the 1 1 U -perturbation theory 
of Harris and Lange [see Eq. (fj-2|)]. Concludingly, the 
IAA improves upon the MAA with respect to the 
weak-coupling and upon the EHA what concerns the 
strong-coupling regime. 



(3) 



VI. DISCUSSION AND RESULTS 
A. Paramagnetic phase 

For the symmetric case of half-filling and paramag- 
netism, electron-hole symmetry requires B a = — 
Tq. This implies that the IAA reduces to the EHA in 
this case and thus cannot contribute to an improved de- 
scription of the Mott transition. We therefore focus on 
the non-symmetric case n ^ 1. Fig. 1 shows the spec- 
tral density of the d — oo Hubbard model on the hyper- 
cubic lattice for U = 4 and n = 0.79 as resulting from 
the MAA, the EHA and the IAA (units are chosen such 
that t* = 1). A non-zero temperature has been taken 
(ft = 1/k-gT = 7.2) to enable a direct comparison with 
the essentially exact QMC result of Jarrell and Pruschkc 
]38| ] which is also shown in the figure. 

For each of the three approximate theories, the 
spectrum mainly consists of the two dominant charge- 
excitation peaks (Hubbard bands). Compared with the 
MAA result, we notice that both peaks are significantly 
broadened within the IAA, their spectral weights be- 
ing almost unchanged. Furthermore, while the MAA 
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Fig. 1. Spectral density of the Hubbard model on the d = oo 
hyper-cubic lattice for on-site Coulomb interaction U — 4, 
band-filling n = 0.79 and inverse temperature /3 = 7.2. Re- 
sults for the modified alloy-analogy (MAA), the Edwards- 
Hertz approach (EHA) and the interpolating alloy-analogy- 
based approximation (IAA) in comparison with quantum 
Monte Carlo data (QMC) from Jarrell and Pruschke pq|. 
All energies are given in units of t* — 1. 



predicts nearly symmetric peak shapes, a clear peak 
asymmetry is observed within the IAA spectrum. The 
broadening effect can partly be traced back to a stronger 
quasi-particle damping in the IAA compared with the 
MAA: The imaginary part of the IAA self-energy is sig- 
nificantly larger at the energetic positions of the Hub- 
bard bands. 

The IAA spectrum is found to be only slightly differ- 
ent from the EHA result: Comparing with the EHA, we 
find the IAA spectral density to be somewhat larger at 
its peak maximum in the lower Hubbard band. There 
is more weight in the lower and less weight in the upper 
peak which is shifted to higher energies compared with 
the EHA. 

Besides the high-energy charge excitation peaks, the 
spectrum obtained by QMC simulations (3^] shows up a 
third structure around the Fermi edge at E — 0, which 
is reminiscent of the Kondo resonance. In the present 
situation, away from half-filling and for finite tempera- 
ture, the resonance is strongly broadened but still clearly 
visible. On the other hand, it is missing in the ap- 
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Fig. 2. The same as Fig. 1, but for smaller band-filling n 
0.57. 



proximate approaches even for T = and increased 
band-filling, although the correct trend can be recog- 
nized in the asymmetry of the IAA (and EHA) lower 
peak. Compared with the QMC result, the MAA pre- 
dicts a too large (quasi) gap between the peaks, while it 
is too small within the IAA and EHA. The IAA (EHA) 
yields a width of the lower peak that is quite close to 
the QMC result. The width of the upper peak, how- 
ever, comes out to large, while it is underestimated by 
the MAA. Passing from the EHA to the IAA, the agree- 
ment with the QMC result is improved very slightly. 

Decreasing the band filling, we find the same quali- 
tative differences between the approximate approaches. 
However, they become less pronounced. This can be 
seen in Fig. 2 where the corresponding results are shown 
for n = 0.57. Comparing with the QMC spectrum, some 
discrepancies also remain for lower n, but the overall 
agreement is certainly improved for all approaches con- 
sidered. This is mainly due to the fact that the Kondo- 
like resonance has disappeared in the QMC spectrum. 

With respect to the EHA the main conceptual 
improvement of the IAA consists in the corrected high- 
energy behavior of the self-energy. This is determined 
by the correlation function B— a . In the strong-coupling 
regime £>_ CT leads to an effective shift of the lower 
Hubbard band n- a B- a [cf. Eq. (|F2])]. This is shown in 
Fig. 3 as a function of n. For the IAA as well as for the 
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Fig. 3. Effective band shift n_ CT _B_ CT for the MAA and the 
IAA and n- a B^\ (Hartree-Fock value, IAA) as a function of 
band-filling n. Hyper-cubic lattice, U = 4, j3 = 7.2, nj = nj . 
Results are symmetric to the n — 1 axis. 

MAA the overall dependence on n is quite similar to 
the results known from the spectral-density approach 
|35] , p9~| as well as from the modified perturbation 
theory [[[9]. We notice that (within the IAA) there is 
a small difference between B_„ and its Hartree-Fock 



value B^l only. Consequently, this implies via Eqs. 
(20) and (^) that for the IAA only minor changes can 
be expected with respect to the EHA. This is consistent 
with the results for the spectral density discussed above. 

B. Ferromagnetism 

Within the self-consistent calculation, the correla- 
tion functions B-„ might become spin-dependent. Ac- 
cording to Eq. ( |l2"| ) this would imply a spin-dependent 
shift of the centers of gravity of the two Hubbard bands. 
Therefore, the band shift B _ G is of exceptional impor- 
tance for the possibility and thermodynamical stability 
of spontaneous ferromagnetic order. 

We consider two different cases where essentially ex- 
act QMC results are available: the hyper-cubic lattice 
with the BDOS given by Eq. @), and the BDOS (|) 
corresponding to an fec-type lattice. A ferromagnetic 
instability is indicated by a divergency of the uniform 
static susceptibility which is calculated within the para- 
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Fig. 4. Inverse static susceptibility \ 1 (inset: x) as a func- 
tion of band-filling n for the hyper-cubic lattice, U = 4 and 
f3 = 7.2. Results for the MAA, EHA and IAA. 



magnetic phase in the limit of an infinitesimally small 
applied field H: x ^ d{n^—ni) / dH\n=a- Jarrell and Pr- 
uschke |38| found a non-diverging susceptibility within 
QMC on the hyper-cubic lattice for various fillings and 
temperatures and interactions up to U — 8. This is 
consistent with the findings of Refs. ]39], |l| where it is 
proved that the Nagaoka state ^d[] is instable for any 
n / 1. Let us mention, however, that a partially po- 
larized ferromagnetic state has been found recently for 
U i ► oo within the non-crossing approximation [|42| . 

The QMC result for U = 4 and (3 = 7.2 shows x to be 
a monotonically increasing function of the filling up to 
n = 1 with a slightly negative curvature, d 2 x/dn 2 < 
]3Sj| . Fig. 4 shows the corresponding results from our 
analytical approaches. First, we notice a zero of x 1 
at n — 0.98 within the IAA. A non-zero magnetization 
is found for 0.98 < n < 1.02 indeed. However, this 
has not been analyzed in more detail since obviously 
(comparing with QMC) it has to be interpreted as an 
artefact of the method. All three approaches considered 
cannot be expected to yield reliable results near to half- 
filling: In this range antiferromagnetic spin fluctuations 
are known to become important which manifest them- 
selves in a narrow Kondo-like resonance at the Fermi 
level. The resonance, however, was shown to be absent 
in the analytical approaches (see above). 



Well below half-filling the IAA predicts a rather fea- 
tureless positive susceptibility which increases with in- 
creasing n and shows up a negative curvature. All this 
agrees with the QMC result. It should be mentioned 
that there are only minor changes in the results for tem- 
perature T h 0. In particular, for the IAA we have 
checked that the artificial ferromagnetic solution is al- 
ways confined to the very small range around half-filling. 
The differences with respect to the EHA result are again 
found to be small except for fillings close to n — 1 where 
the results are unphysical anyway. Finally, we notice 
that within the MAA the functional dependence of x 
on n shows up the wrong curvature for smaller n and 
disagrees in this respect with the QMC as well as with 
the IAA result. 

A non-bipartite lattice as well as a strongly asymmet- 
ric BDOS are known to favor ferromagnetic order. This 
has been demonstrated in a comprehensive way by re- 
cent QMC studies |4^, |44|, [Sand variational approaches 
gfj. The results of UhrigQ and Ulmke @ show a 
ferromagnetic phase to be stable in a fairly extended re- 
gion of the U-n-T phase diagram for the d = oo fec-type 
lattice. Ferromagnetism is favored due to the strong 
asymmetry of the BDOS (|J) with its square-root diver- 
gency at the lower band edge. 

Thermodynamically stable ferromagnetic solutions 
are found within the MAA, the EHA as well as in the 
IAA. Fig. 5 shows the temperature dependence of the 
spontaneous magnetization to = rif — n.j for [7 = 4 
(t* = 1) and n — 0.4. In each case a fully polarized 
state (to = n) is obtained for T = 0. Increasing the 
temperature results in a decrease of m. The magnetiza- 
tion curves are Brillouin-function-like, and the estimated 
critical exponent of the order parameter is (3 = 1/2, i. e. 
the expected mean-field value. The inverse susceptibility 
is also shown in Fig. 5. At the respective Curie tempera- 
ture Tc there is a zero of x _1 as it must be for a second- 
order phase transition. For T > Tc the temperature de- 
pendence of x 1 is almost linear and thereby follows the 
Curie- Weiss law. For the critical exponent of the suscep- 
tibility we obtain the mean-field value 7 = 1. The Curie 
temperature predicted by the IAA (Tc = 0.040) almost 
perfectly agrees with the essentially exact QMC result 
{T c = 0.039 @). 



For larger filling, at n = 0.58, QMC data Uc 



44 



are available for the temperature dependence of to and 



The comparison with the MAA and IAA results is 
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Fig. 5. Spontaneous magnetization m = (n-f) — (nj.) and 
inverse static susceptibility \~ l as functions of temperature 
for the Hubbard model on the d = oo fee-type lattice [cf. 
BDOS in Eq. (§)] for U = 4 and n = 0.4. Results for the 
MAA, EHA and IAA. 



Fig. 6. Temperature dependence of m and x _1 f° r U = 4 and 
n=0.58. Results for the MAA and IAA in comparison with 
the quantum Monto Carlo data (QMC) from Vollhardt et 
al. M3| (diamonds/circles, dashed line: linear fit to x 1 / fit 
with a Brillouin function to m). The thin solid line is a linear 
extrapolation of (the IAA) X~ l {T) to lower temperatures. 



shown in Fig. 6. In both cases a stable ferromagnetic so- 
lution and a second-order transition to the paramagnetic 
phase is found. For high temperatures the inverse static 
susceptibilities obey the Curie- Weiss law. As T i— > Tc, 
however, the IAA inverse susceptibility deviates from 
the linear behavior. This is not observed for lower fillings 
(Fig. 5). The zero of x _1 determines the Curie temper- 
ature Tc = 0.031 while the zero obtained from the ex- 
trapolation of the linear high-temperature trend yields 
the paramagnetic Curie temperature = 0.057. The 
correct (QMC) Curie temperature is found in between: 
T c = 0.05 d|, H|; yet the agreement with respect to T c 
is reasonable. However, the non-linear trend of x 
near Tc is at variance with the QMC result. Also the 
T = magnetization (to = 0.47) comes out too small 
since the Brillouin-function fit to the QMC data for to 
indicates a fully polarized ground state (to = n = 0.58). 

The IAA results for n = 0.58 appear to be less reli- 
able than those for n = 0.4. Still, the IAA turns out to 
be superior compared to the MAA, which yields a con- 
siderably too high Curie temperature (T c = 0.130), and 
also to the EHA: At n = 0.58 the EHA does not yield 
a ferromagnetic solution at all, and we are left with the 
paramagnetic phase only. 



In the discussion of the paramagnet, it has been men- 
tioned that because of the absence of the Kondo-likc 
feature in the IAA spectrum, reliable results cannot be 
expected for fillings close to half-filling. A restricted 
range of validity is also found with respect to the mag- 
netic properties. The IAA is able to predict a reliable 
value for the Curie temperature for fillings up to n ss 0.5 
only. Below n = 0.5, however, the agreement between 
the IAA and QMC is convincing: Fig. 7 shows the cal- 
culated filling dependence of Tc for [7 = 4. The phase 
transition is always of second order. For n < 0.5 the 
IAA result falls into the error bounds of the QMC data 
(except for n = 0.2). Possibly, there is a slight ten- 
dency to overestimate Tc in general. For n > 0.5 the 
Curie temperature steeply decreases with increasing fill- 
ing. Beyond n > 0.61 only the paramagnetic solution 
can be found, while the ferromagnetic region extends up 
to n w 0.88 within QMC. 

The MAA considerably overestimates the possibility 
for ferromagnetism. Over the whole n range considered, 
the Curie temperature is too large by more than a factor 
2. Ferromagnetic solutions exist for < n < 1 (close to 
half-filling we encountered numerical difficulties to sta- 
bilize truly converged solutions). The Tc(n) obtained 
by means of the EHA is similar to the IAA curve. As 
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Fig. 7. Filling dependence of the Curie temperature at U = 
4 within the MA A, EH A and IAA in comparison with the 
QMC result (error bars) by Ulmke Q . 

can be seen in Fig. 7, however, Tq is too small compared 
with the QMC result if n > 0.3, and the ferromagnetic 
region shrinks to n < 0.51 in the EHA, thereby under- 
estimating the possibility for ferromagnetic order. The 
IAA, being a combination of the MAA and the EHA, 
yields a Tc in between and the best agreement with the 
QMC result. 

The interpretation for the different transition tem- 
peratures in the different approaches is the following: 
That there is a too high Tc within the MAA can be 
attributed to the fact that quasi-particle damping is un- 
derestimated compared with the IAA (cf . preceding sub- 
section). Damping effects tend to reduce the density of 
states at the Fermi edge which according to Stoner's 
criterion implies a lower Tc p9[| . The difference in the 
Curie temperatures between the EHA and IAA must 
result from the approximation B_ a t— > B^. The tem- 
perature dependence of the corresponding effective band 
shifts is shown in Fig. 8 for n — 0.4. In the paramagnetic 



phase n^ a B^ a and n_ a B_ a are (almost) temperature 
independent. In the ferromagnetic phase below Tc there 
is a spin splitting of both quantities. As T i— ► the fully 
polarized state is reached, and ni = implies a van- 
ishing effective band shift n^B^ = nj.-Bj 1 ' = 0. On the 
other hand, for T t— > we have Ej (E) i— » which im- 
plies via Eqs. (|) and (§|) that n T B T - n^B^ i-> 0. 
Fig. 8 shows that in both cases, a =|, |, the difference 
n-oB-o — n- a B_^ increases with increasing T and is at 
its maximum for T = Tc, i. e. for the paramagnet. For 
the ferromagnet we can conclude that replacing B— a by 
B^l in the self-consistent solution for the self-energy, 
will result in small changes of the spin-dependent spec- 
tral density only. These are comparable to those dis- 
cussed in the previous subsection for the paramagnet. 
More important, however, within the self-consistency cy- 
cle the effective band shift can have a strong influence 
on the magnetic aspects of the final solution. Therefore, 
using the Hartree-Fock approximation for B- a may be 
insufficient to describe the magnetic properties of the 
system. 

Analyzing the continued-fraction expansion of the 
Green function for T = 0, Uhrig Eol could obtain the 
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Fig. 9. Filling dependence of the magnetization at ksT = 
0.01 and different values of U. The lower critical filling 
n c ,i(U), above which ferromagnetism occurs for fcsT = 0.01, 
decreases with increasing U. The open circles indicate the 
filling un(U) below which the system remains fully polar- 
ized at T = 0: rajv(l.O) =0.36±0.03, njv(1.5) =0.48±0.02, 
un (2.0) = 0.62 ± 0.02 (errors due to the extrapolation to 
T = 0). For comparison the exact results for um by Uhrig 
are also shown (filled circles). 



exact boundary U — Ujsi(n) above which the fully po- 
larized state is stable against a single spin flip for the 
fcc-type lattice. Due to the divergence of the BDOS at 
the lower band edge, the Nagaoka state is locally sta- 
ble even for low fillings n i— > 0, the (in-)stability line 
smoothly ends at Un(0) — 0. To compare with these re- 
sults, Fig. 9 shows the filling-dependence of the magne- 
tization for different U at a low temperature T = 0.01. 
At low fillings the system is paramagnetic. Above a 
critical filling n Cj i({7) a non-zero magnetization is found 
which steeply increases with increasing n and (almost) 



reaches the saturation value m = n. Finally, for higher 
n > n c> 2(U) ferromagnetism breaks down again. The 
higher critical filling n c ^(U) exhibits an unusual U de- 
pendence: It is at its maximum for £7 = 2 and decreases 
for larger U > 2. This is believed to be an artefact of 
the method; we expect the IAA to be able to yield reli- 
able information up to n ~ 0.5 only. On the other hand, 
the U dependence of the lower critical filling is plausible: 
n c ,i(U) monotonically decreases with increasing U. We 
have inspected the IAA result for U — 1 in more detail 
for lower temperatures. Even at T = 0.001 the system is 
not fully polarized for low fillings. We found it difficult 
to access still lower temperatures numerically. Extrapo- 
lating the temperature dependence of the magnetization 
to T = 0, however, we find our results to be consistent 
with a ferromagnetic phase extending to n = + . For 
larger fillings extrapolation to T = yields the value un 
up to which the system remains fully polarized. Fig. 9 
shows that the obtained fillings for U — 1.0, 1.5, 2.0 
are close to the results found by Uhrig pC| |. For the 
larger values of U (U > 2) the IAA is not able repro- 
duce un(U), e. g. a partially polarized ground state for 
U = 4 and n — 0.58 is inconsistent with the results of 
Ref. @ but is found in the IAA (see Fig. 6). 

To estimate to which strength of the interaction U 
the IAA is a reliable approach, we have calculated the 
temperature-dependent magnetization for different U at 
the filling n = 0.5. The results are shown in Fig. 10. For 
each U considered we observe ferromagnetic saturation 
for T > and a second-order transition to the paramag- 
netic phase at the respective Curie temperature Tc{U). 
As a function of U the Curie temperature is found to sat- 
urate with T C (U = oo) w 0.044. We also included the 
predictions of the QMC method for T C {U) in Fig. 10. 
The data for U = 2, 4, 5 are taken from the values of the 
interpolation lines at n — 0.5 in Fig. 2 of Ref. [Q. The 
error bars represent the respective typical error as has 
been estimated from the real data points and errors in 
the vicinity of n = 0.5. There is a convincing agreement 
between IAA and QMC for U = 2. At U = 4 the IAA 
result for Tc still falls into the error bounds, while a 
slightly too low Tc is predicted for U = 5. We conclude 
that for strong interaction U the IAA appears to be rea- 
sonable but less convincing compared with the low and 
intermediate U regime 

The temperature-dependence of the spin-resolved 
spectral densities is shown in Fig. 11. We have taken 
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Fig. 10. Temperature dependence of the magnetization at 
n = 0.5 and different values of U. The QMC data for Tq at 
U = 2,4,5 are taken from the interpolation lines in Fig. 2 of 
Ref. E2]; the bars indicate the typical (QMC) errors. 



n = 0.4 and U = 4 where the IAA should yield reliable 
results. For T = the shape of the | spectral density 
is given by the shape of the free BDOS (§) since the 
magnetization is saturated, m = n, and thus £j = 0. 
The square-root singularity at the lower band edge has 
been slightly smeared out for numerical reasons. We 
have checked that this does not affect the presented re- 
sults for m, x" 1 and Tc- The j spectral density is split 
into the lower and the upper Hubbard band, the en- 
ergetic distance of which is roughly given by U. For 
U = 4 the result of Harris and Lange should apply. The 
energetic shift between the f band and the lower J, Hub- 
bard band (centers of gravity) is read off from Fig. 11 
to be Tii ~ ~ 0.40. On the other hand, the calcu- 
lated effective band shift for the j band is n-j-St = 0.42 
(see also Fig. 8). This is consistent with Eq. ( |l2| ) which 
for the present case and averaged over k simplifies to 
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Fig. 11. Spin-dependent spectral density in the ferromag- 
netic phase for three different temperatures (d = oo fee lat- 
tice, U — 4, n = 0.4). Solid line: majority spin (f). Dashed 
line: minority spin (J.). ksT — 0.043 corresponds to the 
Curie temperature. 



Tii — ^-T — T ■ 

With increasing temperature the upper Hubbard 
band in the f spectral density arises. There is an ex- 
change splitting between the f and J, upper Hubbard 
bands being inverse to the splitting between their lower 
counterparts. For T i— > Tc the spin splitting becomes 
gradually smaller. At the same time we observe a trans- 
fer of spectral weight between the two spin channels that 
speeds up the depolarization of the system. At T = Tc 
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Fig. 12. Spin-dependent spectral density for U — 4 and n = 
0.58. Comparison between the IAA result for fcsT = 0.022 
and the QMC result from Q for k B T = 0.04. In both cases 
the magnetization is m = 0.4. 



both spectra have become identical. Compared with the 
T = case, the singularity at the lower band edge has 
disappeared, and a strong broadening of the spectrum 
can be observed which is due to quasi-particle damping. 

The qualitative behavior of the spectral density as 
a function of temperature, i. e. what concerns the spin- 
dependent centers of gravity, effective widths and spec- 
tral weights of the Hubbard bands, very well agrees with 
the results of Harris and Lange. This implies that the 
higher-order corrections in the t/U expansion that have 
been neglected in Eq. ( |l2| ) can be considered to be small 
(for U = 4). 

At the same interaction strength, but for higher 
hlling there is one QMC result available for the spin- 
dependent spectral density in the ferromagnetic phase 
Q. Unfortunately, the hlling (n = 0.58) and the 
temperature (T = 0.04) where the QMC spectrum has 
been obtained, are not well suited for a one-to-one 
comparison between IAA and QMC since the Curie 
temperature predicted by the IAA is below T = 0.04 
(see Fig. 6). We therefore try a comparison of the 



QMC spectrum with the IAA result at T = 0.022 
where the magnetization is non-zero and equal to the 
magnetization found within QMC (m = 0.4). Fig. 12 
shows the resulting spin-dependent spectral densities. A 
quantitative agreement with the QMC spectrum cannot 
be expected: Due to the temperature dependence of 
the effective band shift n- a B- a , the positions and 
weights of the Hubbard bands cannot be the same as 
those in QMC spectrum at a different temperature. 
Qualitatively, we notice that apart from the Hubbard 
bands there is an additional peak in the f and the J. 
QMC spectrum at£« 0.8 and also a peak at E w —0.4 
in the J. channel which is not seen in the IAA result. At 
zero energy we observe a weak shoulder in the f QMC 
spectrum reminiscent of the Kondo resonance. A (much 
weaker) structure at E — is also found by the IAA. 

VII. CONCLUSION 

Exactly solvable limiting cases represent strong nec- 
essary conditions for any approximate solution of the 
Hubbard model. The central idea of the present paper 
has been to develop an analytical expression for the self- 
energy that is compatible with as much rigorous results 
as possible. From this point of view, the interpolating 
alloy-analogy-based approximation (IAA) can be con- 
sidered as one of the best analytical approaches that are 
available for the Hubbard model in infinite dimensions. 

It avoids the defects of the modified alloy analogy 
(MAA) [^9| which is not able to recover Fermi-liquid 
properties in the weak-coupling regime, and it improves 
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upon the Edwards-Hertz approximation (EHA) |2 
which is at variance with the exact strong-coupling re- 
sults of Harris and Lange. The IAA not only recovers 
the trivial U — and the t = limits but also the respec- 
tive first non-trivial corrections: For small U it is exact 
up to order U 2 and can thus be regarded as a general- 
ization of standard weak-coupling perturbation theory. 
For strong U it yields the correct positions and weights 
of the Hubbard bands which are exactly known from 
perturbation theory in t/U. Formally, the IAA is con- 
structed by a proper combination of the MAA with the 
EHA: The general form of the ansatz for the self-energy 
it taken from the EHA, but following the concept of 
the MAA, the "atomic levels" and "concentrations" are 
determined via comparison with the first four exactly 
known moments of the spectral density. As a conse- 
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quence, the coefficients of the high-energy expansion of 
the self-energy are exact up to order 1/E 2 . Further- 
more, when applied to the Falicov-Kimball model, i. e. 
when suppressing the hopping of one spin species, the 
IAA reduces to the conventional alloy-analogy solution 
which has been shown to be exact in d — oo by Brandt 
and Mielsch. 

The main disadvantage of the method consists in the 
fact that it cannot be motivated from a simple and phys- 
ical concept. Inherently, its justification rests on its in- 
terpolating character only. In particular, the IAA is not 
conserving in the sense of Baym and Kadanoff. This 
must be taken into account when interpreting Fermi- 
surface and related properties. For the symmetric case 
of a paramagnet at half-filling the IAA simply reduces 
to the EHA which predicts the Mott transition but also 
a non-Fermi-liquid metallic phase. Within the EHA and 
similarly for the IAA, the latter turns out to be existing 
in a fairly extended region of the T — phase diagram. 
On the other hand for d = oo there are no indications for 
non-Fermi-liquid behavior in combined studies with it- 
erative perturbation theory and quantum Monte Carlo 
(QMC) simulations §). The non-Fermi-liquid metallic 
phase found in the IAA must be tightly connected with 
the alloy-analogy basis of the approximation and thus 
with the limit of the Falicov-Kimball model since it is 
known that in this limit the (mobile) electrons indeed 
do not form a Fermi liquid Q. 

For an a posteriori test of the quality of the IAA, 
we have compared the results with numerically exact 
quantum Monte Carlo data from Jarrell and Pruschke 
§§ as well as from Ulmke and Vollhardt at al. @ [£§, 
with the results of Uhrig concerning the (in-)stability 
of the Nagaoka state and with the results of the 
MAA and the EHA, mainly for the decisive moderate- 
to strong-coupling regime. 

For the paramagnet on the hyper-cubic lattice, the 
differences between the EHA and IAA results are rather 
small, and judged from the QMC data, the IAA only 
slightly improves upon the MAA and EHA. At fillings 
close to half-filling the IAA fails to reproduce the Kondo- 
like quasi-particle peak in the spectrum, and an artifi- 
cial ferromagnetic solution is found. Reliable results can 
only be expected for smaller n. Here, good overall agree- 
ment with QMC data is found regarding the spectral 
density as well as the static susceptibility. 

Considerable improvement upon the MAA and the 



EHA is achieved for the description of spontaneous ferro- 
magnetic order on a d — oo fec-type lattice. For fillings 
up to quarter filling and for weak and moderate inter- 
action strength U, the IAA is able to predict the mag- 
netic properties and particularly the Curie temperatures 
very reliably. Also the results for strong interactions 
U/t* < 16 at n = 0.5 seem to be still plausible. On the 
other hand, the approximation breaks down for higher 
fillings n > 0.5. The MAA yields a clearly too high 
Tc which is due to an underestimation of quasi-particle 
damping. The EHA is as reliable as the IAA at low 
fillings where it predicts a Curie temperature consistent 
with the QMC result. However, for n > 0.3 the Curie 
temperatures are too low, and a paramagnetic ground 
state is found for n > 0.51 (at U = 4). Thus the EHA 
breaks down at a considerably lower filling compared 
with the IAA. 

In our opinion this is due to the effect of a higher- 
order correlation function B^ a which is of special impor- 
tance for several reasons: Inclusion of B_ a is necessary 
to ensure the correct high-energy behavior of the theory 
as well as the correct moments of the spectral density 
(up to in — 3). This is the main condition to recover 
the exact strong-coupling results of Harris and Lange. 
For U i — ► oo and in the ferromagnetic phase, the cor- 
relation function B- a results in a spin- dependent shift 
n- a B- a of the center of gravity of the lower Hubbard 
band. Thus the stability of ferromagnetic order is inti- 
mately related to the effective band shift. It is known 
that ferromagnetic order cannot be found at all or only 
under extreme circumstances within the Hubbard-I or 
the alloy-analogy solution where T appears in place of 
B-cr in the high-energy expansion coefficients of the self- 
energy. 

The apparent importance of the band shift for 
magnetic order provokes the question whether it is 
possible to construct analytical approaches that on the 
one hand are exact up to the m = 3 moment and on 
the other avoid the deficiencies of the IAA concerning 
its low-energy behavior. For example, up to now there 
is no conserving theory to our knowledge that correctly 
accounts for the m = 3 moment. 
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